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ABSTRACT 


During  the  period  of  this  grant  a  high-speed  photographic,  study. 

has  been  made  of  the  initiation  of  explosive  liquids  (using  sparks 

and  impact)  and  of  single  exystals  of  explosive  solids  (using  shock 

and  hot  wires)*  The  physical  hemogenuity  of  the  material  has  been 

* 

shown  to  have  a  marked  effect  on  its  sensitivity.  Studies  of  the 
fracture  of  inert  solids  have  shown  the  importance  of  specimen  geometry 
and  of  crack  speed  on  the  fragmentation  process.  Both  high-speed 
photography  and  ultrasonic  methods  have  been  used  for  the  fracture 
work.  Ono  of  the  main  advances  has  been  in  the  theoretical  work 
on  thermal  explosion  theory. 
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(a)  Theoretical 

Thermal  Explosion  Theory,  The  effect  of  reactant  consumption  on  critical 
conditions . 

We  outline  below  a  general  theory  describing  the  evolution  of 
temperature  with  extent  of  decomposition  in  an  explosive  mass  suddenly 
introduced  into  a  bjgn  tomporaturo  environment.  Most  of  the  assumptions 
common  to  thermal  explosion  thooxy  are  abandoned  and  oxact  expressions 
for  critical  Frcjik-Kamonotski  (F-K)  temperatures  are  derived.  Our 
method  owes  much  to  that  uf  Adler  &  Enig  (Combustion  and  Flame,  8,  97 > 
1964)  but  is  ^ery  much  more  general. 

The  basic  problem.  A  mass  of  oxplcsivo  of  volume  V,  surface  area  S, 
bar  the  uniform  temperature  T  and  loses  heat  to  its  surroundings  which 
are  a'c  a  temperature  To  at  a  rate  v/hich  is  an  arbitrary  function  of 
(I  -  To).  The  rato  of  heat  generation  duo  to  chomical  reaction  is 
proportional  to  the  product  of  an  arbitrary  function  of  T  and  an  arbitrary 
function  ef  the  fractional  extent  of  decomposition  V'J  •  Whon  doos  tho 
system  exhibit  explosive  behaviour? 

Tho  rato  of  heat  accumulation  is  givon  by: 


* 


(1) 


where  c 7  is  the  density,  t  tho  time,  and  C  tho  specific  hcat/unit  mass* 
The  rate  of  heat  loss  due  to  conductioxi/convectiori/radiation  is 


q2  =  H3(T  -  To)  ^  (T  -  To)  (2) 

where  H  has  tho  dimensions  of  a  hoat  transfer  cooffioiont  and  F^  is 
dimensionloss. 


(1) 


She  rate  of  heat  evolution  duo  tc  chemical  reaction  is 


*  nV  & 

q3  =  "  07(1  at 

where  g  is  the  total  exothermicity/unit  mass* 

The  kinetic  law  is  assumed  to  have  tho  separablo  form 


(3) 


dc 

at  =  “  °°  at 


jr  =  k(T,  c0)  /\  (X) 


to 


whore  c  is  the  concentration  of  reactant  (cQ  its  initial  value)  and 


x-1-»  -i  (5) 

and/\  is  on  arbitiary  function  of  tho  fraction  of  reactant  undecomposed* 
Tho  temperature  dependence  of  tho  initial  reaction  rate  is  assumed 
to  have  the  arbitral^  form 

k(T,  c0)  =  K(o0)  .  F2  (T)  c*p  (J|  )  (6) 

wh.'ro  is  an  arbitrary  function  of  tempo ratu*4 • 

Our  basic  equation  is  tho  law  of  conservation  of  onorgy 


*1  =  *2  *  *3 


(7) 


which  wc  prefer  to  writo  in  dimonsionloss  form  by  making  tho  following 
substitutiens : 

E  tm  RTo  „  .  CTo 


v  = 


RTo 
E 

Vo-  qk(Tot  c0)/o0 


8  =  — 2  (T-To),  ea^,pecki2 
RTo^  h  q 


(8) 


HSeToP^(cTo) 

Equations  (l)  through  (8)  give  our  dimensionless  basic  aquation  tho  form 

BfUl-Y-1  (9) 

A  to 


(2) 


whore  p  is  tho  reciprocal  P-K  reduced  adiabatic  tomporaturo  rise  and  , 
which  no  shall  call  tho  Semenov  criterion,  has  tho  critical  valuo  (2.71828— J"1 
in  tho  special  case  considered  in  Semenov's  classical  theory  for  which:- 
c  =0,  p.O.Asl,  ?!  -P2sl 

Tho  function  (5)  (0),  which  to  call  tho  P-K  function,  is  given  by 


0~®_  0-  0/l+eO  VcTo°)  *2 

*  P2  To(l  +  c8)  (10) 

so  that  (H)  has  tho  simple  form 

®  =  £  e9  +  n-n(l  +  «0  -  (p  -  1)  InQ  (ll) 

when  no  adopt  the  plausible  forms: 

P^  s  const.  X  (T  -  To)p~^;  p  =  const.  (12) 

P2  =  ccnst.  XT®  j  as  const.  (13) 


Tho  treatment  bolow  applios  to  tho  general  form  (10 )  of  tho  P-K  function. 

Porm  (ll)  moroly  illustrates  tho  sort  of  concroto  function  no  havo  in 

mind.  Conventional  thermal  explosion  theory  assumes  form  (11 )  with 

m  =  0,  p  =  1,  e  =  0,  i.o.,  ©  =  9. 

Particular  cases  of  tho  kinetic  function  which  are  of  intorest 
th  n 

aro  arbitrary  n  ordor:  f\  =  X  j  n  ^  0  (14) 

autc catalytic:  /\  =  X(l  -  X)  (15) 

Initial  conditions.  Tho  system  starts  at  somo  arbitrary  tempera¬ 
ture  0Q,  with  no  reactant  consumption 

e  =  e0  at  y  =  o,  (16) 

usually  &0  0. 

A  study  has  boon  made  of  tho  form  of  curvos  (9)  subject  to  (16)  in  tho 


(3) 


)  Plano  and  a  simple  analytical  method  for  finding  close  upper 
and  louor  bounds  cn  the  value  of  V  for  which  explosion  just  fails  to 


occur  has  been  almost  entirely  elaborated  and  will  be  reported  in  detail 
later.  Wo  show  hero  how  to  caiculato  tho  critical  value  of  the  F-K 
temperature .  Curves  (9)  have  an  extremum  when 

-<I) 


A  -v 


-1 


(17) 

and  their  curvature  (by  differentiating  (9)  )  is  given  by 


^»-V-V2A-2o-@z 

d?" 

where  Z  =  (A  V  ~1  (1-8  @d)  +  p°A  x 


(18) 

(19) 


Differentiating  again  T7C  find 

v  p2 ^ "  A'2  °  * 2 -®r  (A-2  o-O) 

whore 


g„1A.y-i  «*-©](» 

QA 


dXj 


and 


af  (A"2  o‘®)  =  ~A_3  0-$  |/N$e  g -  2^ 

The  locus  of  tho  inflexions  of  curvos  (9)  is  given  by 

Z  =  0 

(by  equation  (18)  ).  Or  fr*m  (17) 

A  inf.  =  Aextr.  +  ^Ax  inf.  X(®} 

v/hcro 

"X(®)  =  i  -  Q|H) 


9 


(20) 

(21) 

(22) 

(I9i) 

(23) 

(24) 


(4) 


■<  V  Vr: 


*  V 


Now  usually  0  for  all  X,  so  that  (25)  moons 

frinf .  -  Xoxtr.}  =  ^  <  X  '  (25) 

Tho  ossonce  of  oxplosivo  behaviour  is  that  is  discontinuous, 
i*o.,  1  -  Q(BJg  has  zeros,  so  that  tho  inflexion  locus  has  bronchos 
lying  on  alternate  sides  of  tho  locus  of  oxtroma.  In  all  cases  of 

2 

interest  wo  find  a  region  of  tho  (0,^  )  within  which  both  4r  and  aro 

V  do 

positive,  and  following  Adler  &  Enig  wo  say  that  tho  state  point  must 

enter  this  "explosive  peninsular"  for  an  explosion  to  occur.  Critical 

bohaviour  is  observed  when  a  curvo  (9)  touches  the  inflexion  locus,  i.e., 
2  3 

=  0  for  criticality  (26) 

dh2  dh3 

for  then  tho  state  point  just  fails  to  ontor  tho  explosive  peninsular. 
Sorabining  (26)  and  (18)  through  (22)  we  find 

Y  AAtt.  ,  .1  =  0  (27) 

Ax  (i-e@0)2 

i  ii 

In  gonoral  9  can  be  found  by  eliminating  X  from  (27)  and  (l9i),  but 
cr 

certain  special  forms  of  /\  make  tho  procoss  simpler. 

If  term  I  in  (27)  happens  to  be  indopondont  of  X  thon  (27)  is 
immediately  an  equation  in  0Qr  alone.  (Term  II  however  is  novor  inde¬ 
pendent  of  0  for  forms  of  of  intorost  to  us.) 

If  wo  rogard  A  AxX  =  K  =  const,  as  a  difforontial  oquation  in  A 


we  find  that  the  only  solutions  of  physical  intorcst  arc 


(5) 


A  =  Xa  ;  n  \  O;1 


>  0^—  _  i  =  ri 
*  ’  2  1  "  n 

Ax 


Thus  for  tho  general  n  “  order  reaction  (n  not  nocossarily  intogral), 
and  only  in  the  caso  of  a  kinotic  law  of  this  kind,  wo  find  that  tho 
critical  temperature  is  givon  by 

(1  -  e  (H>0)2  =  n(l  -  e2(3)  ee)  (28) 

and  in  particular  ©  is  independent  of  tho  paramotors  W  ,  p. 

For  tho  standard  functions  (ll)  -  (13)  oquation  (28)  givos  a  quartic 
equation  in  9  which  is  readily  solved  (although  the  resulting  expressions 
aro  cumborsomo). 

The  caso  p  =  n,  m  =  0,  is  especially  simplo  and  yields 

g  -  — .?P.  (29) 

cr  1  -  2pe 

For  most  explosives  we  havo  c  1.  The  solution  of  (28)  and  (ll) 


with  e  =  0  is  simply 


®cr  =  p  "  7  Pn 


=  2p  (2p  -  b)  !  2  jt m  -  n  -  ml  (30)  * 

'  3.c  f  c=0  -■ " * 

Tho  valuos  of  9  givon  by  (28)  no  longor  correspond  to  oxplosion  whon 
cr 

3  -  9  0  q  has  no  soros  (say  becauso  c  is  too  largo). 

Thus  there  is  a  termination  in  explosive  bohaviour  when 

1  -  80e  =  1  +  O20e0  »  0  (31) 

Whon  9  is  oliminatod  from  (31)  wo  thus  havo  an  expression  for  tho  maximum 


Only  valuos  corresponding  to  tho  +  sign  represent  oxplosion 


(6) 


possiblo  voluo  of  c  for  oxplosivo  behaviour, 
give 

c  <  &  hi 2-r-g.).  ffi  ^  o 

in 

.  1 

c  <  J  m  =  0 

Much  stronger  limitations  can  bo  imposed  upon  the  valuo  of  c  by 
porforaing  crude  integrations  of  (9).  Those  will  bo  described  in  tho 
noxt  roport. 

An  important  foaturo  in  our  gonoral  examination  is  the  following. 

Onco  having  located  (®cr>  i^or),  as  wo  ha-vo  by  equations  (l9i)  and  (27), 
wo  can  assign  V  •  (3,  b,  l,  n,  p  and/or  any  other  paramo  tors  that  may 
bo  involved  in  F^  and  and  thon  integrato  backwards  numerically  from 
tho  critical  point  until  =  0  (\j  -  0)  is  roachod.  In  this  way  wo 

detormino 

V  or  ='Vor  (P,  e,  m,  n,  p,  ©0) 
without  resorting  to  the  time-consuming  ’’straddling”  techniques. 

(b)  Liquid  Byplocivos 

Work  undor  this  heading  has  beer,  described  in  tho  papor  ”Formation 
.)f  Cavities  and  Micro  Jots  in  Liquids  and  their  Rolo  in  Initiation  and 
Growth  of  Explosion”  by  Bowden  and  McOnic  (Proc.  Roy.  Soc.  in  Pross). 

In  tho  oxporimonts  a  continuous  access  framing  camera  was  usod 
to  film  tho  oxplosion  process  in  thin  films  of  nitroglycorino.  Initiation 
of  oxplosion  was  by  a  variety  of  methods  including  impact,  spark,  and 
exploding  foil.  Tho  study  showod: 

(i)  that  cavities  or  bubblos  in  the  liquid  could  chango  a  slow  burning 


Equations  (31 )  and  (ll) 


^  for  oxplosion 
) 


(7) 


procoss  into  a  fast  detonation .  Furthor,  that  if  cavities  or  bubbles 
troro  not  present  boforo  tho  procoss  started,  that  thoy  could  bo  intro¬ 
duced  by  stross  waves  sot  up  by  tho  explosion  procoss  itself. 

(ii)  That  thoro  aro  various  situations  during  tho  explosion  procoss  when 
fast  micro- jots  can  bo  produced,  and  that  thesG  last  jots  can  lead  to 
violont  explosion.  Those  jots  oon  bo  producod,  for  oxamplo,  betwoon 
coalescing  drops,  when  cavities  aro  collapsed  by  stress  wavos,  or  whon 
cavities  collapse  near  a  solid  surface. 

Present  work  involves  a  closer  study  of  tho  collapse  of  cavitios  in 
gels  of  ranging  vicsocity. 

(c)  Solid  Explosivos 

Work  on  the  mechanism  of  tho  slow  and  oxplosivo  decomposition  of 
tho  fulminates  was  completed  this  year  and  tho  ro suits  will  bo  published 
shortly.  In  addition  to  normal  methods  usod  to  study  decomposition  and 
measure  oloctron  onorgy  levels  in  solids,  wo  have  carrlod  out  X-ray 
studios  to  dotormino  how  fulminatos  pack  in  tho  unit  coll.  This  was 
required  for  an  understanding  of  the  decomposition  on  tho  molocular  lovol. 

(d)  Fracture  of  Inort  Crystals 

A  study  of  tho  propagation  of  fracturo  in  crystalline  solids  is  in 
progress.  Tho  fracturo  growth  has  boon  followod  both  by  high-spood 
photographic  and  ultrasonic  techniques.  This  lattor  mothod  consists 
ossontially  of  passing  a  high  froquoncy  (say  5  x  10^  o.p.s.)  boom  of 
ultrasonic  wavos  through  a  solid  at  tho  time  of  fracturo.  Tho  osoilla- 
ting  tensions  addod  by  tho  ultrasonic  boom  modulate  tho  fracturo  growth 
oausing  a  faint  ripple  to  bo  producod  on  tho  fracturo  surfaco;  this  acts 


i-:-- 


(8) 


aa  a  porn^aiont  record  of  tho  fracture  dovolopmont,  Sinoo  the  fre¬ 
quency  of  tho  imposod  wavos  is  known  tho  spacing  of  the  ripplos  givos 
an  accurato  mcasuromont  of  crack  velocity,  and  acceleration* 

Tho  fracture  study  has  shown : 

(i)  that  maximum  propagation  ratos  exist  for  many  materials* 

(ii)  that  cracks  in  glassliko  solids  tend  to  branch  whon  they  roach 
a  high  velocity,  hut  that  this  branching  is  suppressed  in  orystallino 
solids  having  well  defined  olcavage  planes,  thu3  enabling  eraoks  in 
oiyctallino  solids  to  roach  relatively  higher  valocitios  (as  compared 
with  their  stress  wave  velocities). 

(iii)  that  '•rack  branching  can  bo  caused  by  tho  interaction  of  a 
fracture  with  a  stress  wave. 

(iv)  that  the  detailed  fragmentation  of  a  solid  is  greatly  affeotod 
by  tho  duration  of  tho  stress  pulsos  involved  in  tho  process. 

Recently  an  attompt  has  been  made  to  photograph  down  through  crystals 
onto  tho  developing  fracture  face.  This  differs  frem  tho  more  conven¬ 
tional  approach  which  observes  crack  propagation  through  solids  edge-on* 
This  now  method  has  boon  shown  to  bo  feasible;  it  is  more  diffioult 
experimentally  but  has  tho  advaritago  that  the  dctailod  development  of 
all  the  features  on  the  fraotvro  faoo  can  bo  followed* 


Roforoncoa 


1*  Adlor  and  Enig,  Combustion  and  Plano,  8,  97,  1964 
2a  Bowdon,  P#  P.  and  MoOnio  M.  P. ,  Proo,  Roy.  Soo.  (in  Pross). 
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